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Abstract 

Refractive processes in strong- field QED are pure quantum processes, which involve only external 
photons and the background electromagnetic field. We show analytically that such processes 
occurring in a plane-wave field are all characterized by a surprisingly modest net exchange of 
laser photons even at ultrarelativistic laser intensities. We obtain this result by a direct calculation 
of the transition matrix element of an arbitrary refractive QED process and accounting exactly for 
the background plane- wave field. A simple physical explanation of this modest net exchange of laser 
photons is provided, based on the fact that the laser field couples with the external photons only 
indirectly through virtual electron-positron pairs. For stronger and stronger laser fields, the pairs 
cover a shorter and shorter distance before they annihilate again, such that the laser can transfer 
to them an energy corresponding to only a few photons. These results apply to both optical and 
x-free electron lasers, and are relevant for upcoming experiments aiming to test strong-field QED 
at present and next-generation facilities. 
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I. INTRODUCTION 



Nonlinear processes have always played a fundamental role in different areas of physics, 
spanning from hydrodynamics, atomic and laser physics to plasma and high-energy physics 
jj. From a theoretical point of view the description of such nonlinear processes, though 
attractive, is also particularly challenging. Since the invention of the laser, it was manifest 
that one of its unique features, the coherence, would allow for the experimental investigation 
of nonlinear phenomena. In a laser beam, in fact, a large number of photons propagate in 
phase and, depending on the laser intensity and on the process at hand, they may act 
cooperatively. One example is atomic high-order harmonic generation (HHG), in which a 
large number of laser photons is absorbed by a single atom and only one high-energy photon 
is emitted (see the reviews 2|)- When laser-driven electrons (mass m and charge e < 0) 
are bound in atoms, nonlinear phenomena start at laser field amplitudes Eq of the order 
of the typical atomic binding field E at = m 2 |e| 5 , which corresponds to a laser intensity of 
J at = E^ t /4iT = 7.0 x 10 16 W/cm 2 (units with h—c—1 are employed throughout). In this 
case the average number of photons absorbed from the laser by the electron is of the order of 
Up/ojQ, where U p = ^E^jmuj^ is its ponderomotive energy and ioq is the central laser photon 
energy. HHG has also been observed for free electrons driven by an intense laser beam, being 
named nonlinear Thomson or nonlinear Compton scattering, depending on if quantum effects 



are negligible or not 



BO 



4j. In both nonlinear Thomson and Compton scattering, the typical 
electric field strength, at which nonlinear effects set on, is given by E re \ = muo/\e\. The 
corresponding intensity is of the order of 10 18 W/cm 2 at optical photon energies ujq ~ 1 eV. 
An electron in a laser field with central laser photon energy uq and electric field strength of 
the order of E Ie \ is accelerated to relativistic velocities already within one laser period and 
its dynamic becomes highly nonlinear with respect to the laser field amplitude [5|. On the 
other hand, quantum effects such as the recoil of the photons emitted by the laser-driven 
electron, strongly modify the emission process when the electric field strength of the laser 
in the initial rest frame of the incoming electron is of the order of the so-called critical 
field E cv = m 2 /\e\ of QED, corresponding to the laser intensity I cr = 4.6 x 10 29 W/cm 2 
6j. Relativistic quantum effects also allow for the nonlinear interaction of a photon with a 
laser field, as in the case of electron-positron pair photo-production (nonlinear Breit- Wheeler 
pair production (NBWPP)) [JQ. This process, as well as any QED process occurring in 



the collision of a photon with a strong laser field (lOj, is essentially controlled by the two 
Lorentz- and gauge-invariant parameters £ = Eo/E ie \ and x = [(kok)/muo]Eo/E CI . Here, 
= ojquj — k ■ k, with k$ = (oj ,k ) and k^ = (u,k) being the four-momentum of 
the laser photons and of the incoming photon, respectively. It is worth observing that in 
the so-called "ultrarelativistic" regime £ 3> 1, the net number of laser photons absorbed in 



for values of £ of the order of 10 2 



NBWPP is very large and of the order of £ 3 Since presently available optical lasers allow 



llj , unprecedented degrees of nonlinearity of the order of 



one million are in principle achievable. 

Refractive QED processes in a strong laser field involve only initial and final photons, and 
the background field [ijj]. Such processes of genuinely quantum nature are a unique tool for 
testing the predictions of strong-field QED on the nonlinear evolution of the electromagnetic 
field in vacuum. Vacuum polarization [13J and photon splitting [14] in a laser field are two 
examples of refractive QED processes, which have been considered in the literature. It has 
been observed in both cases, that the net number of laser photons exchanged with the laser 
field is very small (of the order of unity) even in the ultrarelativistic regime £ ^> 1. 

In the present paper we show analytically that this is a general feature of refractive QED 
processes in a strong laser field. The physical origin of this effect lies in the fact that in a 
refractive QED process, the laser field couples to the external photons only indirectly via 
a virtual electron-positron pair. As we will see below, at higher and higher laser intensi- 
ties the distance covered by the virtual electron and positron before annihilating decreases 
accordingly, in such a way that the process occurs with a net exchange of a low number 
of laser photons. This is in contrast, as we have mentioned, to the NBWPP, which is also 
primed in the collision of a (real) photon and a laser field. However, in NBWPP the final 
electron and positron are on the mass shell, requiring a large amount of laser photons to 
be absorbed for the process to occur at all in the presence of an ultra-relativistic laser field. 
The present results are of relevance for planned experimental campaigns, aiming to measure 
strong-field QED effects in the presence of a background laser field. As we will see, they 
indicate, for example, that, in order to detect nonlinear effects in the laser field amplitude 
on refractive QED processes, it is more convenient to measure the yield of final photons and 
to compare it with the corresponding value without laser field, rather than to measure the 
angular distribution or the energies of the final photons. 
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FIG. 1. Two typical Feynman diagrams relative to a generic refractive QED process in a laser field. 
The thin wavy lines indicate the external photons with four- momenta . . . , and polarization 
four- vectors , . . . , , respectively, and the thick plain lines indicate the laser-dressed electron 
propagators. 



II. CALCULATION OF THE AMPLITUDE OF A GENERIC REFRACTIVE QED 
PROCESS 



Refractive QED processes in a laser field involve in general iVj incoming, N Q outgoing 
photons, with N{ + N Q > 0, and the laser photons. However, for the sake of notational 
simplicity, we consider here the abstract case of only incoming photons (N a = 0) and we set 
Ni = N. The photons have momenta kj and polarization four- vectors , with j — 1, . . . , N 
(see Fig. 1): the jth incoming photon can be "transformed" into an outgoing one via the 
substitutions kj -> — kj and — > * in the amplitude (see Eq. ([T]) below). As it will be 
clear below, the results of the paper are unaffected by this particular choice. The mentioned 
process is described by the sum of all Feynman diagrams, which can be obtained from the 
one in the left side of Fig. 1 by permuting the labels in the photon legs. Among them, we 
consider here only the one in the right part of Fig. 1, and the treatment of the remaining 
diagrams can be performed in an analogous way (any diagram contributing to a refractive 
QED process can always be considered together with the other one, differing only in the 
direction of circulation of the four- momentum through the electron loop). The amplitude 
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M corresponding to the diagrams in Fig. 1 is given by 15] 



M = — e N / d 4 x x ■ • • d 4 x N e- i[(klXl)+ - +{kNXN)] 

J (1) 

x Tv[eiG(xi, x 2 \A)e 2 G(x 2 , x 3 \A) ■ ■ ■ e N G(x N ,x 1 \A)]+ O, 

where the "hat" indicates the contraction of a four-vector with the Dirac gamma matrices 
7 M and where the symbol O indicates the amplitude corresponding to the diagram on the 
right in Fig. 1. In Eq. ([Q) the quantity G(x,y\A) is the dressed electron propagator in the 
laser field. The latter is described by the four- vector potential A 11 = A^((f)), where <fi = (nx), 
with = (1, to) and n being the propagation direction of the laser field. By working in the 
Lorentz gauge, the four- vector potential A 11 ((f)) of the laser field can be chosen in the form 
= (0, A((p)), with n ■ A(<p) = 0. Let ai and a 2 indicate the two possible independent 
laser polarization directions, such that a r ■ a s = 5 rs , with r,s = 1,2, and that a± x a 2 = n. 
Then, the four-vector potential A tl (4>) can be written as A M (</>) = A [aiif>i((j)) + a 2 ip 2 ((p)}, 
where A = —Eq/ojq, aff: = (0, a r ), and the two shape-functions "0 r (0) are arbitrary, smooth 



functions except that they satisfy the relation a/V^ 2 ((f)) + V ; 2 2 (0) — 1 f° r an values of 0, with 
V'i/^^) = c ^V ; i/2(0)/^0- Here, Eq and uo indicate the laser-electric-field amplitude and its 
central angular frequency, respectively [16j . Since the interaction of the jth photon with the 



laser field is controlled by the parameter Xj = r^, with rjj = ujo\kj^x\/'m 2 |6|, EJ, it is natural 
to assume here that ycj ^ for all js, which means k^x 7^ for all js . 

In order to calculate the amplitude M, we employ below the operator technique, devel- 



oped in 



13] for the case of a background plane-wave laser field. In the operator technique 



the electron propagator in the laser field is written as G(x,y\A) = (x\G(A)\y) , where 

G(A) = 1 (2) 
11 — m + te 

with IP = II M (v4) = P^ — eA^((f)) and with e being a positive infinitesimal quantity. Here, the 
four- vector is the four- momentum operator, satisfying the commutation rules P u ] = 
—ig^, where g^ v = diag(+l, — 1, — 1, — 1). By employing the above representation of the 
electron propagator and by using the cyclic property of the trace, the amplitude in Eq. (Tj[|) 
can be simply written as 

M = -e N J d 4 xTr (x^A) ■ ■ ■ G N (A)\x)+ O, (3) 

where we have introduced the block operators Gj(A) = G(A)ij exp[— i(kjx)]. It is convenient 



to express the amplitude M in terms of the "square" propagator 

D(A) = 1 (4) 
IP — m z + %e 

rather than in terms of G(A). The details of the procedure to carry this out are reported in 
the Appendix A. Here, we only provide a summary of this procedure in terms of substitution 
rules. The amplitude M, in fact, turns out to be expressed as 

[N/2]+l 

M = - (M {i) + {1...N ^ N ...1}), (5) 

where are partial amplitudes, with [N/2] indicating the integer part of N/2. The 

quantity {1 . . .N — » iV . . . 1} refers to the fact that each partial amplitude will have N 
indexes corresponding to the N ordered operators Gj(A) in Eq. (j3J), and it indicates that 
the same partial amplitude has to be added, but with the indexes 1, . . . , N appearing 
in the opposite order N, . . . , 1. In turn, each partial amplitude is expressed as a sum 
Ylj'=i Mj of terms and the number Jj of terms in each partial amplitude depends on 
the partial amplitude itself. Each term Mj has the form —e N j d A x Tr (x\Oj^\x), with the 
operator being obtained from the original operator product G\(A) ■ ■ ■ Gjj(A) by means 
of the following substitution rules: 

1. Partial amplitude M^: substitute each block Gj(A) by Dj(A) = D(A) exp[-z(%x)] [2(IIe i )+ 
kjij] (this partial amplitude contains one term). 

2. Partial amplitude combine two successive blocks Gj(A)Gj + \(A) (for j = 
1, . . . , N) and substitute this quantity with the "contraction" 

—Cjj + i(A) = —D(A)ej exp[—i(kjx)]ej + iexp[—i(kj + ix)], then substitute the remain- 
ing blocks as in 1.; it is understood that Gn+i(A) = Gi(A) and that Cn,n+i(A) = 
Cn,i(A); this partial amplitude contains N terms. 

3. Partial amplitude combine twice two successive blocks Gj(A)Gj + i(A) and 
Gf(A)Gf + i(A) (for j = 1, . . . , N - 2, and for f = 3, . . . , N - 1 (if j = 1) or for 
j' = j + 2, . . . , N (if j > 1)), and substitute these quantities with the contractions 
—Cjj + i(A) and — Cyj/ +1 (A), respectively; then substitute the remaining blocks as 
in 1.; it is understood that G^ + \{A) = G\(A) and that Cn,n+i(A) = Cn,i(A); this 
partial amplitude has to be considered only if N > 4 and it contains iV(iV — 3)/2 
terms. 
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4. The above procedure continues by increasing by one the number of combinations 
of successive blocks. The last partial amplitude M^- N ^ +1 > contains the two terms 
(-1)^C 1 , 2 (A)C 3A (A)---C N _ 1>N {A) and (-1) N / 2 C N>1 {A)C 2 , 3 (A) ■ ■ ■ C N - 2 ,n-i with 
N/2 contractions if N is even, or the N terms {-l)^ N ~ 1 ^ 2 D 1 {A)C %3 (A) ■ ■ • C N -i tN (A), 
(-1)(^-^C NA (A)D 2 (A)C 3 4A) ■ ■ - C^^A),..., 

(-l)( N -V/ 2 C lj2 (A)C 3A (A) ■ ■ ■ C N ^ N ^(A)D N (A) with (N - l)/2 contractions if N is 
odd. 



in 



Now, a useful exponential representation of the square propagator D(A) has been found 
see also 



D(A) = — i [°° dse is{tl2 ~ m2+ie) = -i I™ dse- i{m2 - it)s \l + ^[A(<p + 2sP x ) - M<f))]\ 
Jo Jo I 2P X J 



x e -*/o ds'[P±-eA{4>+2s' P X )\ 2 ^isP^Px ^ 



(6) 



where we have introduced the operators P^ = (P t + P x S)/2 and P x = — (Pt — Pxn) = —{nP) 
of the conjugated momenta to the coordinates = t—x\\ and X = (t+x»)/2, with x» = n x, 
such that and X can be interpreted as a "time" and a "space" coordinate, respectively, 
i.e., {<p, P+] = -% and [X, P x \ = i. Note that t = X + <p/2, x\\ = X - <p/2, P t = P^ - P x /2, 
and P H =P^ + Px/2. 

Out of the different partial amplitudes which arise from the above substitutions, we work 
out only the following one 

M« = -e N J d 4 xTi (x\D(A)e~ l{klx) [2(U ei ) + Mi] ■ ■ ■ D(A)e- l{kNx) [2(Ue N ) + k N e N )\x), 

(7) 

which arises from the substitution in 1.. This partial amplitude is always present, inde- 
pendently of the number of the external photons and, as it will also be clear from the 
considerations below, the analysis of the other partial amplitudes proceeds analogously. By 
looking at the expression of the operators D(A) (see Eq. (jHJ)), the coordinate operators X 
and x± appear to occur only in the exponentials relative to the external photons. By em- 
ploying the operator identity e t<ykjxS> f (P)e~ t<ykjX ' > = f(P + kj), we can move all the operators 



e i(k j , x x+k jz ±-x±) ^ Q left and let them act on the bra (x\. The result is 

M« = - e N J d 4 xe i{KxX+K ^ x ^Tr (x\e~ ik ^ {2[(W + ^)e liM ] + k 1 e 1 }D 2 (A) 

• • • x e - lk »-^{2i(W + K%)e N ^J + k N ^e N ^}D N (A) ( 8 ) 
x e- ikN ^[2(Ue N ) + k N e N ]D(A)\x), 

where K» = £f =1 tf, = £jV fcf (note that < = if"), and D t (A) = D{A)\ Px ^ Px+Kl x ^ P±+Kl ^ 
with I = 2, . . . , N. Now, the operators between the bra (x\ and the ket \x) do not contain 
the coordinates X and xj_, and the identities 

(X\f(Px)\X) = J ^f(Px), (x ± \g(P±)\*±) = J J^9(P±) (9) 

valid for arbitrary functions f(Px) and g(Pj_) can be applied (we assumed here that the 
eigenstates \p) of the four-momentum operator P M , i.e., P M |p) = p^\p), are such that (x\p) = 
e -i(px) an j (p\p'} = {2n) A 5 i {p — p')). Moreover, the integrals in X and x± are easily taken 
and the partial amplitude becomes 



M (1) = - (-ie) N 5(K x )5 2 (K ± ) j d<p j dp x J d 2 p ± J d Sl ■ ■ ■ ds N e 
x Tr (0|{2[(^ - eA((P))e N J + k N e N } 
x 1 1 + ^n[A(<f> + 2s lPx )) - A(<j>)} | 

x {2[(p^-eA(0) + <) ei J + Mi} 
x { 1 + 777 ^ rn[i(0 + 2s 2 (p x + k 2 , x )) - A(0)] 1 

X e~ 1 ^ 2 ds 2[P-L+ K 2,±-eA(</>+2s^(px+K2,x))] 2 e -2is2P0(px+K2,x) e -«2,^ 

• • • x {2[(p" - eA(0) + «w_i)ejv-ij + ^v-ie^-i} 



— i(m 2 — ie)S 



(10) 



x ^ 1 + T< T W + 2sjv ^* + " 

I 2(Px + KjV,x) 

X e~ l ^ )N ds 'NiP-L+ K N,±-eM<P+ 2s ' N (Px+^N,x))] 2 e -2is JV P^ ) (p x +K ]Vi x) e -«JV,00|0^ 

where 5 = Si+- • -+sat. We note that in this expression of the amplitude, we have substituted 
the operator P M with the number p^ + k 1 * in the four- dimensional scalar products (Pej). 
First, we observe that, since (kjej) = 0, then it is (ftjej) = (Kj + iej), for j = 1, ... ,N — 1 
and (knCn) = 0. Moreover, although the substitution (Pej) — > (p^ + Kj)ej tfi is evident for 

8 



the components px and p± (see Eq. flU]) and the definition of the operators Di(A) below Eq. 
(jHJ) ) , it is in principle not justified for the remaining component P^. However, we show in the 
Appendix B that gauge invariance implies that the four- dimensional scalar products (Pe 3 -) 
actually do not involve the component P^. The remaining matrix element can be calculated 
by employing the identity 

e-^|0) = |0-0o), (11) 

where 0o is a constant, and the fact that {4>\4>') = S((p — <fi'). The resulting ^-function 
5(2si(px + Kijc) H — " + 2sn(px + Kn,x)) can be exploited to perform the integral in px and 
the result is 

M« = - 5(K X )6 2 (K ± ) Jd^J d 2 P± J™ dSl ' - dSN e -**-*)s e -<K** 

x Tr /{2[(p" - eA(4>))e N J + k N e N } (l + — —h[A(<f> + 0i) - A((f>)} 



25n h 



x 



x {2[(p" - eA(<f) + 00 + Kl) ei J + Mi} \ 1 + 1 r^—n[A((f> + 2 ) - A(<j> + &)] 
• • ■ x {2[(p M - eA(0 + ^.x) + ^_ 1 )e JV -i, At ] + kN-ie N -i} 



x U + — n[A(0 + <p N ) - A(cj) + 

(12) 

In this expression we have simplified the notation by introducing the "average" 

i N r s * 

/=cE ds 'ifM) (13) 

of iV arbitrary functions fj(s'j), the residuals 



and the quantities 



and 



0j = 2 ^ (15) 

i=l 

^{^s'^rf-eA^ + ^l (16) 
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with 



b[ = 26k 1<x s[ (17) 



1-1 



# = 2 ^ ^i,xsi + 25k 1jX s' u l = 2,...,N. (18) 

i=l 

Note also that px = —Rx, that 4>n = and that in our gauge irj t x/<j>{<j>i s'j) — K j,x/4>- 
Moreover, in Eq. ( IT2|) and in the successive expression of M^ 1 ' , the quantity px in the trace 
has to be interpreted as —Rx- 

In order to take the integral in p±, it is convenient first to shift p± as p± — > p± — 
7T±(0> { s }) ; where {s} = Si, . . . , s^v- In this way, the resulting expression of the amplitude 
can be written as 

I N 

x Tr( IJ{2[(p^ + Stftf, Sj )) ej J + kjej} (19) 
\ j= i 

x { 1 + 2^^ "^ + ~ i(0 + ^ } ) ' 



where 



n<t>, {s}) = J2 / ds 'M^ s 'j) 5 ^m # 2 

,=1 ^0 



m 2 + ze], (20) 



where Sk n+ i = 5k\ and 4>x+i = </>i- The integral in p^ = (pi,P2) can be written as a sum 
of integrals of the form 

I ni ,n 2 = J d 2 p xP r P T e~ iSp2 \ (21) 

where n\ and n 2 are two non-negative integers. The integral J ni)Tla vanishes if ni and/or n 2 
are odd, whereas it is equal to 

(m-1)!!^-!)!! 

Jni ' na ^ (2i5-)(ni+n2+2)/2 

if ni and n 2 are both even. In conclusion, we can write the partial amplitude in the 
compact form 

M (i) =tl(-ie) N 6(Kx)6 2 (K ± ) J d<j) J°° e -i[K^- F{M s})] 
I N 

x Tr( l[{2[(p^ + 5^(0, Sj ))e^\ + %e,} (23) 



x |l + 2 ^ 6 +i x n[A{4> + ~ ^(0 + 0i)] |y 
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where the substitution rules 



(24) 

(poi) \ 11 / (pa 2 ) JO ^ n i an d/or n 2 are odd 




. . . . (25) 

V-°?/ VV-«iy 1 ^ggfef if m and n 2 are even 

in the expression of the trace are understood. 

Before passing to the estimation of the net number of laser photon exchanged in a refrac- 
tive QED process, we observe here that the integral representation 

-i poo 

T]- 2 \ = (-W d Sl ---ds N e^^ ds >'- m2 ^\ (26) 

~~ Pj -m 2 + ie J 

of the electron propagator in vacuum in momentum space, suggests to interpret the quantity 
5iTj((j),s'j) as an "effective" instantaneous four-momentum of the virtual particle flowing 
between the (j — l)th and the jth vertex (see Eqs. (|23|) and ([20]) ). 



III. ESTIMATION OF THE NET NUMBER OF EXCHANGED LASER PHO- 
TONS 

If there were no external laser field, the remaining integral in in Eq. (|23|) would provide 
the 5-function 8{K^) : which, together with the other three 5-functions, would imply the 
overall energy-momentum conservation = 0, as expected. In the presence of the laser 
field, a measure of the net number of photons exchanged with the laser field during the 
refractive QED process is determined by the quantity K^/uo, where ujq is the central laser 
angular frequency. In order to estimate the net number of laser photons exchanged, we recall 
that the multiphoton nature of the process is controlled by the parameter £ = \e\E /mui , 
where Eq is the amplitude of the electric field of the laser 6|, |7| . From the physical meaning 
of this parameter, in fact, it is not surprising that if £ < 1, the net number of photons 
exchanged with the laser field is of the order of unity. Thus, we directly consider below the 
ultra-relativistic case where £ 3> 1. In order to further specify the physical regime, we have 
also to consider the parameters (see the discussion below Eq. ([[])). If largely exceeds 
unity, an electron-positron pair can be in principle created in the collision of the laser field 
and the jth external photon. The subsequent emission of radiation by such a pair would 
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represent a background for the refractive QED process. Thus, we limit here to the case where 
the parameters Xj are fixed and less or of the order of unity, such that electron-positron pair 
production from laser-external photons is negligible. Correspondingly, we also exclude the 
possibility that electron-positron pairs can be created only by the external photons, even 
though, as it will be clear below, the following considerations will not depend formally on 
this condition. 

It is convenient to write explicitly 

6rf(<f>, fySn^te, s'j) = -26K j)X 6KM - [Stt^, s'^f (27) 

and to shift the variable (f) as <fi — > <fi + $, with $ such that 

JV 

+ 2 SK 3,xSKj,<t> 8 3 = °- ( 28 ) 

In this way, the the partial amplitude can be written in the convenient form 
M « =^(- ie ) N 5(K x )5 2 (K L ) J d<f> y°° ds ^- s2 dsN e-W+F^As})] 
I N 

x Tr/ H{2[(pP + 5^(0 + $, Sj ))e^\ + kjej} (29) 
\ j= i 

x j 1 + r-^ n[A(4> + $ + <f> j+1 ) - A(<f> + $ + <p 3 )} 

where 




The advantage of this form with respect to that in Eq. ( 1231) is that all the N integrands in 
+ {s}) are strictly positive and therefore that F^((j) + $, {s}) > 0. This implies, in 
fact, that the integration region in dsi ■ ■ ■ ds^ mainly contributing to the partial amplitude 
is confined to sufficiently small values of Sj such that that F^((f) + $, {s}) < 1, as 
otherwise the function exp(— iFj,{<p + {s})) would be highly oscillating. From what we 
mentioned at the beginning of this section, this would already indicate that the net number 
of photon exchanged during the refractive QED process is of the order of unity. However, 
in order to complete the proof, we have still to analyze the pre-exponential function. In 
fact, if N is small, then the different powers of the external field present in this function 
would not essentially change the net number of laser photons exchanged. However, this 

12 




could in principle occur for large Ns. In order to show that this is not the case, we recall 
that in the considered regime, the parameters rjj = are much smaller than unity and 
therefore, in the effective integration region with respect to the variables si, . . . , sjvj h is 
^o|^i,x| s j ^ Uo\5nj^x\/'m 2 "C 1, where we used the fact that Sj < 1/m 2 (see Eq. (I3"U|) ). 
Consequently, it results that w |<Aj|, ^ol^l ^ 1 an d> by assuming that j/c^l < \K^\ for all 
js, that ujo\&\ *C 1 (see Eq. (128]) ). This observation allows one to expand the four- vector 
potential in Eq. f )29|) as [l8| 

+ $ + 0j) _ 2^(0) ($ + <K**) (31) 

^ i=i ' 

^(0 + $ + 0;.) - 2^(0) ( $ + + *«i.* s i) > ( 32 ) 

where £^(0) = — tL4 M (0)/d0 (note that £^(0) is not a four- vector) . Analogously, one obtains 

i(0 + $ + i+1 ) - i(0 + $ + 0,-) « - 2E(<j))5K j+ljX s j+1 (33) 

i-i 



5^(0 + $, 4) + 2e^(0) 



b 1=1 \ <=1 2 / 



(34) 



Now, the fact that £^(0 + $, {s}) < 1 implies, as an order-of-magnitude estimate, that 
[^ 7r j,±(0 + s j)] 2s j ~ 1/-/V. Thus, the above expansions, together with the fact that ~ 
1/yS (see Eq. ( 125]) ). indicate that in the effective formation region of the process, the ratio 
between the terms in the pre-exponent proportional to the laser field and those which do 
not contain the laser field itself is less than unity. Therefore, terms containing higher powers 
of the external field are subdominant and, in conclusion, the probability of an exchange of 
a net number of photons much larger than unity is suppressed also for large values of N. 

In order to make our analysis more concrete, we consider the particular case of a 
monochromatic, circularly polarized laser field. In this case, the vector potential is given by 
A(0) = — (Eo/cu )[cos(cj (f>)ai + sin(co> o 0)a 2 ]. Starting again from the general expression in 
Eq. f )29|) (see also Eq. ( 130]) ). it is convenient to introduce the vectors 

a iA s 'j) = C j( s j) a i + S j ( y s' j )a 2 (35) 
a i)S (4) = -$(4)0! + 0,(4)02, (36) 
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where Cj(s'j) = cos(o;o(^ ) + 4>'j)) an d Sj( s 'j) = sin(wo( < ^ ) + <j>'j))- 111 this way, we obtain 

<$7Tj,_L(0 + s'j) = 5nj ; ± - m£[cos(u} (fi)Saj !C (s'j) + sin(a> o 0)<fab>(4)] (37) 
and the function F^cj) + $, {s}) can be written as 

F^{<f) + $, {s}) = F {{s}) + F c {{s}) cos(uW>) + F s {{s}) sm{u <j)), (38) 

where 

Fo({s})=J2 / ^.((^, x ) 2 + m 2 {l + e 2 [(^(4)) 2 + (^(^)) 2 ]} - (39) 

N 



^c/ s ({s}) = - 2m£ / ds'jSKj^ ■ 5a j:C/s 

.7=1 J ° 



>;•)• ( 4 o) 



Note that the integrals in ds'j in Fo({s}) and F c / S ({s}) can be easily taken in the present 
case, which is however not necessary here. The discussion below Eq. (13"U1) indicates that in 
the effective integration region it is F ({s}), |-F c /s({s})| < 1. We consider now the prototype 
integral in 

1({ S }) = J d0 e -^ +i ^ + *' {s})1 , (41) 

which is present in the partial amplitude M' 1 '. After introducing the quantities Fa({s}) 
and ifo({s}) according to the definitions 

F c ({s}) = F A ({s})cos(M{s})), (42) 
F.({s}) = F A (W)Bin(W)(W)). ( 43 ) 

and after passing to the variable <p = u (j) — ip ({s}), we obtain 

oo 

l({s}) = 27 re- i [( x ^°^« s » +Fo « s »] i' ni S(K^- ni u )J ni (F A ({s})), (44) 

ni=— oo 

where we employed the identity exp(iz cos (p) = ^^ = ^^i n Jn{z) exp{in(p) in terms of the 
ordinary Bessel functions J n (z) of integer order n, valid for an arbitrary complex number 
z [ji3]. Equation (14"4"1) shows that indicates the net number of photons absorbed from (if 
rii < 0) or ceded to (if n\ > 0) the laser field. The well-known property of ordinary Bessel 
functions J n (x) of a real (positive) argument of being much smaller than unity at n > x 



and the fact that Fa({s}) = yFf(Js}) + F^({s}) < 1 shows, at least for the terms in the 
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pre-exponent independent of the laser field, that the net number of photons exchanged with 
the laser field is of the order of unity. The general observation below Eq. ( I34p indicates that 
also high-order terms in the laser field in the pre-exponential will not essentially increase 
the net number of laser photons exchanged during the refractive QED process. 



IV. DISCUSSION 



As we have already mentioned above, it is interesting to compare the low net exchange 
of laser photon in a refractive QED process with what happens in the case of the NBWPP, 
which does also occur in the collision of a (real) photon and a laser field. In particular, we 
limit again to the strong-field limit corresponding to £ ^> 1 at fixed invariant parameters 
Xj ~ 1. The real electron and positron created via the NBWPP at £ ^> 1 are already 
ultra-relativistic and a large net number of laser photons of the order of £ 3 are absorbed 
from the laser field in order to fulfill energy- momentum conservation [7| . On the other hand, 
a refractive QED process occurs via a virtual electron-positron pair and this manifests itself 
in the appearance of the integrals in dsi ■ • ■ dspf in the partial amplitude At larger and 

larger values of the electric field amplitude, the effective integration region in ds\ ■ ■ -ds^ 
reduces accordingly, in such a way that the function + $, {s}) is always of the order 
of or less than unity, and then that the net number of laser photons exchanged is of the 
order of unity, too. More specifically, we recall that if p M = (e,p) is the momentum of a 
classical electron at the initial value <fi = (A(0) = 0), then the component p</>(0) of the 
four-momentum = (e((j)),p((j))) at is given by [5j] 

U s m 2 + [p ± - eA(0)] 2 

?w) = 5- • ( 45 ) 

By performing the change of variable <fij = 2dKj t xs'j in Eq. fl30l) . we see that F^cfi + $, {s}) 
qualitatively corresponds to the quantity YljLiJ^ 3 1 ^ ( t ) 'jPj,<t>{4 > 'j)i where O — and where 
(<//■) is the component (f) of the four-momentum of the virtual electron/positron flowing 
between the (j — l)th vertex and the jth vertex. Thus, the condition F^{(p + $, {s}) < 1 
corresponds to the fact that, according to Heisenberg uncertainty principle, the virtual 
electron-positron pair annihilates after an interval A<fi'j in <pj given by ~ l/V^j, where 
Vfaj indicates the order of magnitude of the momentum flowing between the (j — l)th 
vertex and the jth vertex. This corroborates the interpretation that in a refractive QED 
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process, the stronger is the laser field, the higher is the four-momentum flowing through 
the electron-positron loop. Accordingly, the virtual electron-positron pair propagates for a 
shorter distance inside the laser field, such that the net number of photons, that can be 
exchanged in the process is always of the order of unity. 

This difference between the net number of photons exchanged with the laser field in a 
general refractive QED process and in NBWPP could appear at first sight to contradict 
the optical theorem, when the imaginary part of the amplitude of a refractive QED process 
can be related to the total rate of the corresponding pair-production process (e.g., the 
refractive QED process corresponding to NBWPP is essentially the polarization operator) 



15] . However, this is not the case, because the total rate of a pair-production process 
does not contain information on the net number of photons exchanged with the laser field, 
as the rate is integrated over the whole phase space of the created electron and positron. 
More quantitatively, since a plane-wave field depends only on the spacetime variable 0, it 
is possible to write th ^-matrix element Sfi of an arbitrary process occurring in such a 
background field as 

S ft = 5 fl + z(27r) 3 5 2 (P /i± - Pi,A_)5{P fiX - P i>x )R f i, (46) 



where indicates the total initial/final four- momentum. The optical theorem 15] here 
reads 

2Im(i2«) = ^(2tt)V(P /i± - P,±)5(P f , x - Pi,x)\Rfi\ 2 (47) 
/ 

and we are interested to the case in which in the initial state there are a certain number 
of photons, whereas in the final state an electron-positron pair is present. By limiting, for 
simplicity, to the case of a monochromatic laser field with angular frequency Uo, we can 
expand the amplitude Rfi as 

oo 

R fi = ( 27r W P /.* - p i,i> - niUo)T niJi , (48) 

n;=— oo 

and the optical theorem provides the relation 

oo 

2Im(T 0i44 )= ^(2vr) 4 5(P;-Pr-^o^)|T niJl | 2 . (49) 

ni=-oo f 

On the one hand, this identity shows that only the quantity T 0iii corresponding to no net 
exchange of laser photons in a refractive QED process is relevant for the optical theorem. 
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On the other hand, as already mentioned, all the quantities \T ni ji\ 2 corresponding to the net 
exchange of an arbitrary number of laser photons in the pair-production process are to be 
included on the right-hand side of Eq. (I48p . In the specific example of NBWPP, the above 
conclusion is confirmed by the fact that the total pair production rate at £ ^> 1 becomes in- 
dependent of the parameter £ (it depends only on the parameter x = (u>o\kx\/m 2 )(E /E cr ), 
where k^ is the four- momentum of the external photon), and it coincides with the corre- 
sponding total rate in a constant-crossed field but averaged over the laser phase [7|. 

It is also worth observing that, although the net number of laser photons exchanged in a 
refractive QED process is of the order of unity, high-order terms in the laser field amplitude 
contribute to the process (as, for example, in the Bessel functions in Eq. ()44p ). Such 
nonlinear terms stem for the exchange of laser photons without a net absorption or emission 
during the process. The fact that + $, {s}) < 1 (that F A ({s}) < 1 in Eq. (jUJ) for 
the case of a circularly-polarized, monochromatic laser field), suggests that in general the 
exchange of a large number of laser photons is not probable. However, such nonlinear effects 
can strongly modify the amplitude of a refractive QED process. This observation suggests 
that, in general, in order to detect nonlinear effects in the laser amplitude in a refractive 
QED process, it is more convenient to measure yields of final photons and to compare them 
with the expected result without laser field, rather than to measure, for example, the energy 
or the angular distribution of the final photons. In fact, the optimal regime of parameters 
to detect a refractive QED process is at Xj ~ 1. Now, even considering next generation 
of 10-PW optical laser systems [6j], providing an intensity of the order of 10 23 W/cm 2 , the 
ratio Eq/E cx is smaller that 5 x 10~ 4 . Thus, in order to have Xj ~ 1, initial photon energies 
are required of the order of 1 GeV. For final photon energies of this order of magnitude, 
if only a few photons from an optical laser (ujq ~ 1 eV) are effectively exchanged, it is 
not feasible in practice to measure the induced effects on the final photons' energies and/or 
angular distribution. The same conclusion can also be drawn in the case of present and 



under-construction x-ray laser facilities [6], for which even the ultra- relativistic limit £ 3> 1 
is not feasible, due to the relatively large photon energy (w > 1 KeV). 

The analysis carried out so far accounts only for the tree-level diagrams like those in Fig. 
1, i.e., it has been assumed that radiative corrections are negligible. In the presence of an 
external plane-wave field this is the case if olx 2 -^ z <C 1 for all j, where a = e 2 ~ 1/137 is the 
fine-structure constant, i.e., if Xj <C 1/a 3 / 2 « 10 3 [7J]. However, radiative corrections and 
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high-order diagrams would in any case involve only virtual particles in such a way that the 
physical argument given above and concerning the net number of laser photon exchanged 
would again apply On the other hand, as we have already mentioned, the regime Xj 3> 1 is 
not suitable for observing a refractive QED process, due to the background photons emitted 
by the produced electron-positron pairs. 



V. CONCLUSION 



In conclusion, by employing the operator technique, we have shown in general that refrac- 
tive QED processes in a laser field occur with a net absorption/emission of only a few laser 
photons even in the ultrarelativistic regime £ 3> 1. Due to the general validity of the above 
analysis, the present investigation is of relevance for experimental campaigns at upcoming 
optical and x-ray laser facilities. Our main conclusion is that in order to experimentally 
observe nonlinear effects in the laser field amplitude in such processes, it is more convenient 
to measure yields of final photons in a refractive QED process and compare them with the 
expected value without laser field, rather than, for example, to measure the energies or the 
angular distribution of the final photons. 
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APPENDIX A 



In the present appendix we will indicate how to express the amplitude ([3]) in such a 
way that it contains only the square propagators D(A) (see Eq. (jl])). It is convenient 
to introduce here the notation (note that some of the above symbols have been already 
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introduced between Eq. ([3]) and Eq. (JB])) 



Cjj + i(A) =D(A)ij exp[—i(kjx)}ij + i exp[—i(kj+ix)]. 



Gj(A) =G(A)ej exp[-i(kjx)], 

D 3 (A) =D(A) e X p[-i(%x)][2(n ei ) + kjej], 

Qj(A) =D(A)ej eotp[-i(A; i a;)]Gf- 1 (A), 



(50) 
(51) 
(52) 
(53) 



The following identities can be easily proven 



G S (A) =D j (A)-Q j (A), 
Qj{A)Dj+i{A) =Qj{A)Q j+ i(A)+C jtj+1 (A), 



(54) 
(55) 



where for j = N, the index N + 1 has to be intended as 1 (recall the cyclic property of the 
trace). In order to further simplify the notation, we also define the generalized trace of a 
matrix operator O 



Since, as will be clear, the procedure to transform the quantity Tn(A) only depends on if 
N is odd or even, we explicitly work out only the cases N = 3 and iV = 4, being the cases 
N > 4 completely analogous. Now, 



T 3 (A) = Tr x [G 1 (A)G 2 {A)G 3 (A)}+ O 

= Tr x [(Di(A) - Q 1 (A))(D 2 (A) - Q 2 (A)){D 3 (A) - Q 3 (A)))+ O 

= Tr x [D 1 (A)D 2 (A)D 3 (A)]-Tr x [Q 1 (A)D 2 (A)D 3 (A)] 

- Tr x [D 1 (A)Q 2 (A)D 3 (A)} - Tr x [D 1 (A)D 2 (A)Q 3 (A)} + Tr x [D 1 (A)Q 2 (A)Q 3 (A)} 

+ Tr x [Q 1 (A)D 2 (A)Q 3 (A)\ + Tr X [Q x ( A) Q 2 (A)D 3 ( A)] - Tr x [Q 1 (A)Q 2 (A)Q 3 (A)}+ O . 



The first term already contains only square propagators and, by applying the identity (1551) 
to the three terms containing only one operator Qj(A), we see that the contributions com- 
ing from the first term in Eq. fl55|) exactly cancel the terms containing two operators 




(56) 



such that it is sufficient to analyze the quantity 



T N (A) = Tr x [G 1 {A)-..G N (A)]+0 . 



(57) 



(58) 
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Qj(A)Qj+i(A). Thus, we obtain 



T 3 (A) = Tr x [D 1 (A)D 2 (A)D 3 (A)] - Tr x [C 1>2 {A)D 3 (A)} - Tr x [D 1 (A)C 2 ^(A)] 

(59) 

- Tz x {C^{A)D 2 {A)\ - Tr x [Q l (A)Q 2 (A)Q 3 (A)}+ . 
Now, we consider separately the quantity 

T +i3 (A) = Tr x [Q 1 (A)Q 2 (A)Q 3 (A)] = Tr x [G + (A)e 1 e- j(fe ^ ) G + (A)e 2 e^ (fc ^ ) G + (A)e 3 e- j(fc ^ ) ] , 

(60) 

where we have introduced the quantity G + (A) = (II + m + ie)^ 1 , which corresponds to the 
electron propagator but with m — > —m. By imagining to work in the Dirac representation 
of the gamma matrices 15j, we consider the unitary matrix U = j ^ 2 ^ 5 and we note 
that U^^U^ = 7 M '*, where the upper index t indicates the transpose with respect to the 
Dirac-matrices indexes. Since the four-momentum operator is hermitian, it is easy to show 
that UG + (A)W = — [G(— A)] tx , where the upper index t x indicates the transpose with 
respect to the Dirac-matrices and to the spacetime indexes. In this way, by inserting the 
unity operator UU' in Eq. ( 160]) before and after each ij and by exploiting the fact that 
Tr^O^Of) = Tr x [(0 2 Oi)' a: ] = Tr x (0 2 Oi) for arbitrary operators 1 and 2 , we obtain 

T +>3 (A) =Tr a [Q 1 (A)Q 2 {A)Q 3 (A)] = -Tr x [G(-A)e 3 e- J ( fc3 ^G(-A)e 2 e^( fc2;c )G'(-A)e 1 e- i ( fcl ^] 
= - Tr x [G 3 {-A)G 2 (-A)G 1 (-A)}. 

(61) 

Now, we recall that, in general, the quantity Tn(A) also contain the contribution from the 
Feynman diagram where the electron arrows are reversed (see Fig. 1) and that, due to Furry 
theorem [1^], only terms proportional to an odd power of laser amplitude effectively con- 
tribute to T 3 (A), i.e., T 3 (A) = — T 3 (— A). Therefore, by applying the same above procedure 
to the additional contribution from the Feynman diagram where the electron arrows are 
reversed, we obtain 

T 3 (A) =^-{Tr x [D 1 (A)D 2 (A)D s (A)] - Tr X [C 12 (A)D 3 (A)] 

2 (62) 

- Tr^AfaM] - Tr x [C 3jl (A)D 2 (A)} + {123 -> 321})}, 

where the quantity {123 — > 321} means that the previous terms have to be added, but 
with the indexes 1, 2 and 3 appearing in the opposite order 3, 2 and 1. This result exactly 
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corresponds to the general procedure given in the main text below Eq. fl3]) for the case 
N = 3. 

The case with N = 4 can be worked out in a completely analogous way and we only 
stress the differences with respect to the case N = 3. The starting point is the quantity 

T 4 (A) = Tr x [G 1 (A)G 3 {A)G 3 (A)G i (A)]+ O 

= Tr^D^A) - Q 1 {A))(D 2 {A) - Q 2 {A)){D 3 (A) - Q 3 (A)){D 4 (A) - Q 4 (A))} + O 
= Tr x [D 1 (A)D 2 (A)D 3 (A)D 4 (A)} - Tr x [Q 1 (A)D 2 (A)D 3 (A)D i (A)) 

- Tr x [D 1 (A)Q 2 (A)D 3 (A)D 4 (A)} - Tr x [D 1 (A)D 2 (A)Q 3 (A)D 4 (A)] 

- Tr x [D 1 (A)D 2 (A)D 3 (A)Q 4 (A)\ + Tt x [Q 1 (A)Q 2 (A)D 3 (A)D 4 (A)} 
+ Tr x [Q l (A)D 2 (A)Q 3 (A)D 4 (A)\ + Tr x [Q 1 (A)D 2 (A)D 3 (A)Q 4 (A)\ 
+ Tr x [D x {A)Q 2 {A)Q 3 {A)D 4 {A)) + Tr x [D x {A)Q 2 {A)D 3 {A)Q 4 {A)) 
+ Tr x [D x {A)D 2 {A)Q 3 {A)Q 4 {A)) - Tr x [D x (A) Q 2 (A) Q 3 (A) Q 4 (A)] 

- Tr x [Q l {A)D 2 {A)Q 3 {A)Q 4 {A)} - Tr x [Q 1 (A)Q 2 (A)D 3 (A)Q 4 (A)} 

- Tr x [Q 1 {A)Q 2 {A)Q 3 {A)D 4 {A)} + Tr x [Q 1 (A)Q 2 (A)Q 3 (A)Q 4 (A)}+ . 

(63) 

By applying the identity fl53|) in the terms containing only one operator Qj(A), four of the 
six terms with two operators Qj(A) and Qji(A) cancel, and we obtain 

T 4 (A) = Tr x [D 1 (A)D 2 (A)D 3 (A)D 4 (A)] - Tr x [C 1;2 (A)D 3 {A)D 4 (A)\ 

- Tr x [Di (A) C 2) 3 (A) D 4 (A) ] - Tr^ [D x (A) D 2 (A) C 3)4 (A) ] - Tr x [C 4il (A)D 2 (A)D 3 (A)} 
+ Ti x [Q l {A)D 2 {A)Q 3 {A)D 4 {A)\ + Ti x [D l {A)Q 2 {A)D 3 {A)Q 4 {A)} 

- Tr x [D 1 (A)Q 2 (A)Q 3 {A)Q 4 (A)} - Tr x [Q 1 (A)D 2 (A)Q 3 (A)Q 4 (A)] 

- Tr x [Q 1 (A)Q 2 (A)D 3 (A)Q 4 (A)} - Tr x [Q 1 {A)Q 2 (A)Q 3 (A)D 4 (A)] 
+ Tr x [Q 1 (A)Q 2 (A)Q 3 (A)Q 4 (A)]+ o . 

(64) 



By applying the identity (1551) in the remaining terms containing two operators Qj{A) and 
Qj'(A), two of the four terms with three operators Qj(A), Qj>(A) and Qj»(A) cancel, and 
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we obtain 



T 4 (A) = Tr x [D 1 {A)D 2 (A)D 3 (A)D 4 (A)) - Tt x [C 1 , 2 (A)D 3 (A)D A (A)] 

- Tr x [D l {A)C 2 , 3 {A)D A {A)\ - Ti x [D 1 (A)D 2 (A)C 3A (A)} - Tr x [C Atl {A)D 2 {A)D 3 {A)\ 
+ Tv x [C 1 , 2 (A)Q 3 (A)D 4 (A)} + li^D^C^Q^A)} 

- Ti x [Q 1 (A)D 2 (A)Q 3 (A)Q 4 (A)} - Tr x [Q 1 (A)Q 2 (A)D 3 (A)Q 4 (A)} 
+ Ti x [Q 1 (A)Q2(A)Q 3 (A)Q,(A)]+ . 

(65) 

Finally, by applying again the identity (155]) in the two terms containing two operators Qj{A) 
and Qf(A), the new terms containing three operators Qj(A), Qj>(A) and Qj»{A) combine 
to the remaining two terms also containing three operators Qj(A), Qji(A) and Qjn(A), and 
give two terms Ty x [Q 1 (A)Q 2 (A)Q 3 (A)Q i (A)] with a minus sign. In conclusion, we have 

T 4 (A) = Tr x [D 1 (A)D 2 (A)D 3 (A)D 4 (A)) - Tt x [C x , 2 {A)D 3 {A)D A {A)} 

- Tr x [D x {A)C^A)D^A)\ - Ti x [D 1 (A)D 2 (A)C 3A (A)} - Tr x [C Atl {A)D 2 {A)D 3 (A)\ 
+ Tr x [C 1;2 {A)C 3A {A)\ + Tr x [C^{A)C 2t3 {A)\ - Tr x [Q 1 (A)Q 2 {A)Q 3 (A)Q i {A)}+ O . 

(66) 

The trace Tr x [Qi(A)Q 2 (A)Q 3 (A)Q4 i (A)] can be manipulated exactly as in the case N = 3 
and we arrive to the final result 

T 4 (A) = ^{Tr x [D x {A)D 2 {A)D 3 {A)D A {A)\ - Tr x [C 1>2 {A)D 3 (A)D 4 {A)\ 

- Tr x [D l {A)C 2 , 3 {A)D A {A)] - Ti x [D 1 (A)D 2 (A)C 3A (A)} - Tr x [C Atl {A)D 2 {A)D 3 {A)\ 
+ Tt x [C 1>2 {A)C 3A {A)] + Tr x [C 4>1 (A)C 2)3 (A)} + {1234 -> 4321}}, 

(67) 

which again corresponds to the substitution rules given below Eq. ([3]) for the case N = 4. 
APPENDIX B 

In this appendix, we show that the four- dimensional scalar products {Pej) do not contain 
the operator P^. We temporarily assume that kj 2 ^ for all js. In this way, by introducing 
the quantities fj? u = n^a" — n u a^ with r = 1, 2, the four-vector can be expanded with 
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respect to the basis [1 



j ~ — FT - ' J ' ~ — FT - ' ( ^ 



A (3),M = A (4)„ = (6g) 

as ej = Eti^A^, with fef = -(Af ej ) (note that (Af A< v) ) = -5 UV , with u,-y = 
1, . . . , 4). If we write the total amplitude M as M = ej )W • • • ejv,^M M1 "' MAr , then 

M= ^ ^•••feS^^-.A^M^™ (70) 

ni,...,iijv=l 

and gauge invariance requires that k\^ 1 M IJ,1 '"^ N = ■ ■ ■ = k^^ N M IJ,1 '"^ N = [l5|. This implies 

(3) 

that the terms proportional to the four- vectors A^- ^ . and those proportional to the divergent 
part of the four-vectors in the limits kj 2 — > 0, do not contribute to M. Thus, the 
amplitude M remains finite in the same limits kj 2 — > 0. Moreover, the quantities (Pe/) only 
effectively involve contractions of P M either with n^ 1 or with a±, 2 , so that they do not contain 
the operator P^. It is also worth pointing out here that in the limit kj 2 — > 0, although 
the contributing part of A^ goes to zero as yj kj 2 , the corresponding contribution to the 
amplitude M remains finite because the quantity tip = — (A^e 3 -) diverges as 1/ J k 2 in the 
same limit (see Eq. fl69|) ). In conclusion, by means of the above limiting procedure, our 
analysis can also be applied to the case in which any of the external photons is real, i.e., 
on-shell. 
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